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Abstract 

We derive the microscopic spectral density of the Dirac operator in SU{Nc > 3) Yang- 
Mills theory coupled to Nf fermions in the fundamental representation. An essential 
technical ingredient is an exact rewriting of this density in terms of integrations over 
the super Riemannian manifold Gl{Nf + 1|1). The result agrees exactly with earlier 
calculations based on Random Matrix Theory. 
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1. Introduction 


For the understanding of chiral symmetry breaking (and restoration) in QCD and 
other QCD-like theories with massless or near-massless quarks it is essential to know the 
distribution of the small eigenvalues of the Dirac operator. In the conventional thermo¬ 
dynamic limit where the space-time volume V is taken to infinity prior to taking quark 
masses to zero, the accumulation of eigenvalues A near A ~ 0 may or may not lead to a 
non-vanishing density of the eigenvalues at the origin, p(0)/D. According to the Banks- 
Casher relation |l[], Sq = = '^p(0)/V, this spectral density is an order parameter for 

chiral symmetry breaking. In all that follows we only consider the case where the chiral 
condensate does not vanish in the chiral limit. 

The full spectral density p(A) is prohibitively difficult to compute, sensitive as it is to 
physics at all scales from the infrared up to the ultraviolet cut-off. Because of its relation 
to the chiral condensate, even just computing p(0) is tantamount to understanding in all 
detail the dynamics that underlies chiral symmetry breaking in QCD. Moreover, since 
p(0) involves an energy scale in the deepest infrared, only non-perturbative techniques 
such as lattice regularizations can have a hope of computing this single number from 
first principles. The fact that only low-momentum modes of the Dirac spectrum are of 
relevance here nevertheless opens up the interesting possibility of studying the problem 
in a new setting, using the low-momentum representation of the QCD partition func¬ 
tion. This effective theory is a chiral Lagrangian. Its equivalence with the conventional 
representation of QCD in terms of microscopic degrees of freedom (quarks and gluons) 
becomes exact in precisely the limit of interest here: zero (or almost-zero) momentum^. 
There is only one input parameter in this approach: the chiral condensate Sq, and hence 
the value of spectral density evaluated at the origin, p(0). For computational purposes it 
is convenient to impose a condition which ensures that only exact zero-momentum modes 
dominate the euclidean QCD partition function. This can be achieved by restricting the 
(large) space-time volume V to obey the inequalities 

^QCD 

where m.^ is the mass of the Goldstone bosons associated with the lightest quark mass 
1^, Q. The hrst inequality is required to separate the Goldstone modes from the other 
hadronic excitations with a mass scale of Aqqq or higher. When D —> cx) the second 
inequality gives an unphysical limit in which the pion is ultra-light, and never fits into the 
space-time volume. However, a wealth of information about the Dirac operator spectrum 
can be computed exactly in this limit. 

^The domain of validity of the zero momentum mode approximation is determined by the pion decay 
constant F |^. 


( 1 ) 
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Once it is assured that only zero-momentum modes contribute to the effective QCD 
partition function, its evaluation becomes remarkably simple. The space-time integration 
of the effective Lagrangian density yields just an overall volume factor, and the partition 
function becomes identical to a zero-dimensional group integral over the coset determined 
by the pattern of chiral symmetry breaking. As was noted by Leutwyler and Smilga |^, 
it is highly advantageous to consider these partition functions in sectors corresponding 
to definite topological charge z/ (in order to avoid absolute value signs, we will take v 
positive or zero from now on). Consider SU{Nc > 3) gauge theories with Nf fermions 
in the fundamental representation. If chiral symmetry breaks according to SUiiNf) x 
SU]i{Nf) —>• SUv{Nf), the effective partition function of the zero momentum modes takes 
on the form ||, ^ 

dU det^{U)exp[VJ:oRe Tt[MU^]} (2) 

in a sector of topological charge u, and with fermion mass matrix Ad. Note that the 
effective partition function only depends on masses rrii and four-volume V in the scaling 
combination of /ij = rUjldSo. A few years ago it was realized that this effective partition 
function itself has an entirely different representation in terms of large-A^ Random Matrix 
Theory [^, The precise relationship is as follows. 

Define a “chiral” random matrix partition function with the same global symmetries 
as the field theory partition function by B0 

Nf 

/=i 


where 


V 


0 ikh \ 
iW^ 0 j ’ 


( 4 ) 


and W is a nx m matrix with u = |n — m| and N = n + m. The Random Matrix Theory 
potential V{W^W) is not restricted by any symmetry conditions, but it can be shown that 
all relevant results in the proper limit do not depend on this potential once one imposes 
that the matrix spectral density at the origin p(0) be non-vanishing B a 0 o, 0, 

The equivalent of the topological charge u in the effective partition function (|^) is taken 
to be fixed. In the large-A^ limit we thus have n = N/2. The matrix elements of W 
are either real (/3 = 1, chiral Gaussian Orthogonal Ensemble (chGOE)), complex (/3 = 2, 
chiral Gaussian Unitary Ensemble (chGUE)), or quaternion real {j3 = 4, chiral Gaussian 
Symplectic Ensemble (chGSE)) [0. For Yang-Mills theory with three or more colors and 
quarks in the fundamental representation the matrix elements of the Dirac operator are 
complex, and we have (3 = 2. It can be demonstrated that in the microscopic domain (|l|) 
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the random matrix partition function for any potential V{W^W) can be mapped exactly 
onto the effective hnite volume partition function (Q). This was shown in ref. ^ for a 
Gaussian potential, and, because the random matrix partition function does not depend 
on ViW'^W) in this limit [jT^, it holds for any potential. For more discussion of the 


random matrix partition function we refer to |15 


Because of the spontaneous breaking of chiral symmetry, the smallest eigenvalues of 
the Dirac operator are spaced as l/p(0) = vr/SoD. In order to study the behavior of the 
smallest eigenvalues in the approach to the thermodynamic limit, it is natural to rescale 
the eigenvalues according to u = AI/Sq and to introduce the microscopic limit of the 
spectral density 


1 


u 


Ps{u) = lim 

V^oo P Ln FLn 


(5) 


For broken chiral symmetry this results in a nontrivial limiting function. 

To understand the signihcance of the identihcation between the chiral Lagrangian 
(Q) and the Random Matrix Theory (|]) in the microscopic domain, it helps to view 
these partition functions as generating functions for the chiral condensates {'ipi'ipi). 

Conventionally they are dehned by taking the zero-mass limit in the end. However, it is 
useful to focus instead on the mass-dependent chiral condensate, dehned in the obvious 
way by 


• • • 7 PNf') — 


1 d 


log Z^[ni, . . . , /iTVf] • 


( 6 ) 


As is evident from the spectral representation of the condensate, this quantity carries 
much information about the microscopic spectral density of the Dirac operator: 


^ . s ^ , Ps{u; fxi,..., fXNf) 

LR/ii,... ,/rjvJ = ^0 du -^^^ 

■' J lU + fXi 


(7) 


The issue at hand is whether this relation can be uniquely inverted to provide the mi¬ 
croscopic spectral density ps{u] ni,..., unj) in terms of the mass-dependent chiral con¬ 
densate Sj(pi,...,/ijv^). While the relation (0) resembles the Stieltjes transform of 
Ps{u; pi,..., fiNf) (which under suitable convergence criteria has a unique inverse), the 
/Xj-dependence of the microscopic spectral density ruins this identihcation. This problem 
was recently solved in a paper by three of us [|^, where it was noted that the introduction 
of an additional fermion species into the theory can be used to provide the needed unique 
inverse of a relation of the form (|^). Clearly what is needed is that the microscopic spec¬ 
tral density becomes insensitive to the addition of this additional species. The solution is 
to simultaneously introduce yet another quark species, but this time of opposite statistics 
[T6[| . In the original held theory formulation this corresponds to a Euclidean partition 
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function of the form 





det(ip — m^i) 
det(ip — my2) 


det(i^ — ruf) 
/=i 


^-Sym[A] 


( 8 ) 


When rriyi = m ^2 this partition function simply coincides with the original one. However, 
it is now also the generator of a mass-dependent chiral condensate for the additional (say, 
fermionic) quark species, i.e.. 


• • • ,mAr ) = 


d 


V drriyi 


logZP^ 


(9) 


— lTly2 — 


In terms of the spectral density the valence quark mass dependence of the chiral conden¬ 
sate can be rewritten as 


• • ■ ,mNf) 


= i— 

V k + ruy 


1 

V 


d\ 


p(A;mi, • • ■,mNf) 
iX -|- rriy 


( 10 ) 


Here, (• ■ ■) denotes an average with respect to the distribution of the eigenvalues. No¬ 
tice that the spectral density in this equation is for rriyi = my 2 and coincides with 
the original QCD one. Contrary to (I), the relation (p!0D can then be inverted to give 
p(A;mi,... ,mjv^). As mentioned in |^, the spectral density follows from the disconti¬ 
nuity across the imaginary axis. 


Disc|„^^.;^S(m^) = lirnE(iA + e)-E(iA-e) = 27r ^(5(A + A^)) = 27rp(A), (11) 

k 

where we have suppressed the dependence on the sea-quark masses. 

The formulation of the spectral density as a derivative of the partition function (^) 
provides a natural explanation for the result that the microscopic spectral density can be 
related to the usual hnite volume partition function with two additional flavors p0| , pTf . 

There is a close analogy to the quenching prescription of lattice gauge theory, but 
it is worthwhile stressing that, since we eventually restrict ourselves to equal masses 
rriyi = m„ 2 , there are no approximations involved in introducing two additional quark 
species of opposite statistics in this manner. 

The effective Lagrangian corresponding to the partition function (|D with two addi¬ 
tional quark species is determined by its underlying global (super-)symmetry structure. 
Its precise form is very close to the usual one, except for the fact that the group manifold 
of the Goldstone modes is that of a supergroup. The main object of this paper is the 
calculation of the valence quark mass dependence of the chiral condensate from this effec¬ 
tive theory. Below we will focus on the domain (Q) where the nonzero momentum modes 
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factorize from the effective partition function. This super-symmetric effective Lagrangian 


has been considered in a different context by Bernard and Golterman 1^^ who, by 
analogy with lattice gauge theory, call it the partially quenched chiral Lagrangian. The 
terms “sea quarks” for the original physical fermions and “valence quarks” for the ad¬ 
ditional (fermionic) species have been borrowed from lattice gauge theory as well. In 
perturbation theory the cancellation of the associated determinants in the path integral 
corresponds to the omission of loops with this fermion. We stress again that in the present 
context this is not an approximation, but precisely what is required for the associated 
spectral density to be equal to the QCD spectral density. The quenched version of this 


effective Lagrangian was recently derived from a two-sub lattice random flux model 
using the flavor-color transformation introduced by Zirnbauer 


Below we will find that the microscopic spectral density obtained from the extreme 
infrared sector of (H) agrees with the result BH derived from chiral Random Matrix 
Theory. One can wonder how exact results can be obtained this way. The answer is 
simple. In the microscopic scaling regime (P the effective QCD Lagrangian is entirely 
free of dynamics. The only assumption is that the theory supports spontaneous breaking 
of chiral symmetry. With this knowledge alone the microscopic limit of the QCD partition 
function can be written in terms of a zero-dimensional group integral over the Coldstone 
manifold. This is both true for the usual chiral Lagrangian and for the chiral Lagrangian 
corresponding to a partition function with valence quarks. The latter results in the exact 
distribution of the eigenvalues of the Dirac operator in the microscopic scaling limit. This 
computable and universal end of the Dirac operator spectrum is, as was to be expected, 
entirely independent of the detailed QCD dynamics. 


The organization of this paper is as follows. In section 2 we introduce the effective 
theory corresponding to QCD with additional quark species of different statistics and 
discuss the geometry of the Coldstone manifold. Some of the pitfalls associated with the 
evaluation of super-integrals are discussed in section 3. In section 4 we calculate the mea¬ 
sure using an explicit representation of the super-unitary group. In section 5 we evaluate 
the valence quark mass dependence of the chiral condensate in the sector of topological 
charge u, and derive directly from this the microscopic spectral density Ps(A), which is 
found to agree exactly with earlier results based on Random Matrix Theory. A general 
expression for arbitrary number of flavors and topological charge is obtained in section 
6. This expression is evaluated in section 7 for an arbitrary number of massless quarks 
in the sector of zero topological charge and for one sea quark in a sector of topological 
charge u. Concluding remarks are made in section 8. Additional technical details of the 
calculations can be found in appendices A and B. 
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2. Low Energy limit of QCD 


If chiral invariance is spontaneously broken in QCD, the low energy limit of the theory 
is dominated by the Goldstone modes associated with this spontaneous symmetry break¬ 
ing. The basic flavor symmetry of the QCD action with Nf (physical) sea quarks and 
valence quarks is 


GhiNf + iV.liVQ ® GlniNf + N,\N,). (12) 

The reason for this bigger symmetry group is that a priori we do not relate the integration 
variables in the partition function by complex conjugation. This symmetry group is not 
necessarily a symmetry of the QCD partition function. The symmetry transformations 
may violate the convergence of the integrals. There are no problems for the Grassmann 
integration. However, the integrations of the bosonic quark fields are only convergent if 
the helds are related by complex conjugation. For this reason a Ua{Nv) transformation 
on the bosonic quark helds is not a symmetry of the partition function. However, a 
Gl{Ny)/U{Ny) axial transformation is consistent with convergence requirements. 

After spontaneous breaking of the chiral symmetry according to 


ShiNf + Ny\Ny) 0 SlniNf + Ny\Ny) ^ SlviNf + Ny\Ny), 


(13) 


the symmetry of the QCD partition function is reduced to Slv{Nf + Ny\Ny) 0 GZy(l) 
where 0 denotes the semi-direct product. Notice that an axial Gl(l) subgroup of the 
group (|l^) is broken explicitly by the anomaly. 


The Goldstone manifold corresponding to the symmetry breaking pattern (|^ is based 
on the symmetric superspace SIaINj + Nv\Ny). In our effective partition function the 
terms that break the axial symmetry will be included explicitly resulting in an inte¬ 
gration manifold given by SlA^Nf + Ny\Ny) 0 G/^(l). However, this manifold is not a 
super-Riemannian manifold and is not suitable as an integration domain for the low en¬ 
ergy partition function. As an integration domain we choose a maximum Riemannian 
submanifold of G/^(iVj -|- Ny\Ny). This results in a fermion-fermion block given by the 
compact domain SUa^Nj -f Ny), whereas the boson-boson block is restricted to the non¬ 
compact domain Gl{Ny)/U{Ny). Because of the super-trace, this compact/non-compact 
structure is required for obtaining a positive definite quadratic form for the mass term 
and the kinetic term of our low energy effective partition function [^. For a detailed 
mathematical discussion of this construction we refer to a paper by Zirnbauer m- 

In this paper we restrict ourselves to the case of just one valence quark, Ny = 1, which 
is all that is needed to derive the microscopic spectral density of the QCD Dirac operator. 
We will denote our integration manifold by Gl{Nf -|- 1|1). The fields are parametrized by 


U = exp(iA/2$/F), 


(14) 
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where $ is a {Nf + 2) x {Nf + 2) superfield: 


$ 



0 

X 



( 15 ) 


and F is the pion decay constant. Here, 0 is a {Nf + 1) x [Nf + 1) Hermitean matrix 
containing the ordinary mesons made of quarks and antiquarks. The factor i in the field 0 
provides us with a parametrization of GI{1)/U{1). It represents a meson made out of two 
ghost quarks. Finally, y and y (not related by complex conjugation) represent fermionic 
mesons consisting of a ghost-quark and an ordinary anti-quark. They are the Goldstone 
fermions associated with the spontaneous breaking of the chiral supersymmetry (^). We 
choose a diagonal mass matrix with Nf sea quark masses resulting in the (A^/-|-2) x {Nf+2) 
quark-mass matrix 


A4 = diag(mi,..., my — J). 


(16) 


The source J acts as to lift the degeneracy between the valence quark and its superpartner. 
We will eventually set J = 0. 

The axial symmetry of the QCD partition function in the sector of topological charge v 
is broken explicitly by a mismatch in the number of left-handed and right-handed modes. 
Under a global axial transformation Ua, the integrand in (|]) is multiplied by a factor 
SdeG{UA)- Taking into account the explicit breaking by the mass term as well, the QCD 
partition function in the range O reduces to the effective low energy partition function 

n 

(A^) = / dUSdeU(U)e^^s*'’(^^+^^“0. (17) 

^ JueGl{Nf+l\l) 

As discussed before, the integration manifold is given by the maximum Riemannian 
submanifold for the symmetric superspace Gl{Nf + 1|1). For a definition of the super¬ 
determinant, Sdet, and the super-trace, Str, we refer to the book by Efetov |p8| . 

A consistency check of this effective partition follows by comparing ( 0 ) with (^, 
which implies 

l.,=o = Z«h(V<) ■ ( 18 ) 

While this identity may appear quite trivial, only a careful analysis of the superintegral 
involved on the left hand side of this equation will ensure that it is fulfilled. This issue is 
connected with the appearance of Efetov-Wegner terms, which we will discuss next. 


3. A Simple Example of a Superintegral 





In this section we remind the reader of some of the problems that occur in integra¬ 
tions over a supermanifold. In spite of the fact that Grassmann integrations are always 
convergent, the actual integration over a supermanifold can be quite an involved task. 
Even an apparently innocent change of variables may result in a subtle paradox |^, 


To illustrate this fact, let us look at an example given by Zirnbauer 0. Consider the 
supermatrix 


A = 


a a 
(3 ib 


and the following Gaussian integral 

X = j da db da djS = j da db da df3 . 


(19) 


( 20 ) 


By explicit integration or by Wegner’s theorem |^, one knows that 1=1. After Zirn¬ 
bauer, let us change variables according to 


A' = 


s — (Tr(s — it) —r(s — it) 
cr(s — it) it — aT(s — it) 


The Berezinian of this transformation is easily computed to be: 


d{a, b, a, j3) 


Therefore the integral 


a(s,t,cr, r) 
becomes 


it)2' 


( 21 ) 


( 22 ) 


X = ds dt da dr - -r— = ds dt da dr - 


ity 


ity 


(23) 


The change of variables (|2T|) has removed all the Grassmann variables from the integrand. 
One might naively conclude that X = 0, in flagrant contradiction with the result we 
derived before with the original set of variables! However, this conclusion is wrong: The 
Berezinian is singular at s = t = 0, and the integral over s and t in (^31) is therefore not 
dehned. 

An even more puzzling substitution is to keep the fermionic variables a , f3 oi (|^, 
while using the eigenvalues of A as bosonic integration variables |^. They are easily 
computed to be: 

aP 

(24) 


a = a + 
ib = ib + 


a—ib 

Oi0 

a—ib * 


The Berezinian of this change of coordinates is equal to one, and the integral ( PP]) therefore 
becomes 


X = da db da dp e 2 


-i(a2+fe2) 


(25) 
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Again the Grassmann variables have disappeared from the integrand, and one could 
naively conclude that X = 0. This is not true. The substitution above generates boundary 
contributions, the so-called Efetov-Wegner terms pSj |^ . In order to compensate for the 
nilpotent terms in the integration domain of the superintegral, the measure has to be 
modihed accordingly. The correct measure, including the Efetov-Wegner term, is known 
in this example . 


This short discussion exemplihes the subtleties involved in the computation of a super¬ 
integral. In a superintegral, one has to be suspicious about any change of variables. When 
bosonic coordinates are shifted by nilpotent terms, the path of integration may contain 
even functions of the Grassmann variables. This can result in the notorious Efetov-Wegner 
terms in the measure pM |3^. This is not an academic problem. In 


our case, if the integration of the Gl{Nf -|- 1|1) partition function is performed via a su¬ 
persymmetric generalization of the Itzykson-Zuber integral |5^, pp|, Efetov-Wegner 
terms appear in the measure in a way very similar to the second change of variables in 
the example above. However, as is evident from the example, the possible appearance 
of Efetov-Wegner terms depends on the choice of parametrization. Below we introduce a 
parametrization without such anomalous terms for the observables under consideration. 
For example we will End that the partition function is properly normalized without the 
need for anomalous terms. 


4. Parametrization of Gl{Nf + 1|1) and calculation of the measure 


In order to construct an explicit parametrization for Gl{Nf -l- 1|1) we remind the 
reader of a parametrization of the Riemannian superspace associated with G(l|l) which 


appeared earlier in the literature [pS 

/ gie 
0 


0 


exp 



(26) 


The usefulness of this parametrization lies in the factorization into an ordinary and a 
Grassmannian factor. As was discussed in previous section the Goldstone manifold is 
a Riemannian superspace that requires a parametrization in terms of compact and non¬ 
compact variables P0| , p^, 0. More specihcally, instead of the parametrization in 


boson-boson block is given by GI{1)/U{1), obtained by the replacement expi(j) 
The generalization to Gl{Nf + 1|1) is immediate: 


U = 



0 \ 

/ Un 

(o- 

g. ) exp 


0 cxi \ 


Pn 


CXfi 

0 y 


, the 
exp s. 


(27) 
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where 0 is a null-vector of length n = Nf + 1, Un is a unitary matrix and the {«*} and 
{/5j} are Grassmann variables. If this product is written as 

U = UMg, (28) 


we hnd that 


SU' = U-^dUU-^Uo = U-^dUo + dUM-^. 


T-l 


- 1 , 


T-l 


(29) 


Since the Jacobian of a similarity transformation is unity, the invariant measure of the 
super-Riemannian manifold is given by 


d[u\ = B(U„ U,)DU„ds n dPkdat, 


(30) 


k=l 


where DUn is the invariant measure of U{Nf + 1), and the Berezinian is given by 

The differentials SUn are defined by SUn = U~^dUn- From the dehnition of the su¬ 
perdeterminant we hnd that the Berezinian factorizes as follows 


B2 


with 


Bi = Sdet 


5U' 




(32) 


(33) 


and 


Bo = det 


5Ua, 1 ■ ■ ■ JGo 

yln+1 9nn-\-l yn+l 1 yn+ln 

Sai■■■SanSPi ■■ ■ SPn 


(34) 


Here, SUg = dUgU^. By inspection one hnds that the matrix in Bi contains a block 
of zeros. One then trivially verihes that Ri = 1. This implies the factorization of the 
measure in an ordinary unitary invariant measure and a Grassmannian factor. 

The calculation of the second determinant proceeds as follows. We hrst rewrite Ug as 






X 


(35) 


where G denotes the exponent in Ug = exp G and 


X 


Nf+l 

^ = X) Ptat. 


(36) 


k=l 
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Here and below, sir\hx/x and coshx are meant as a formal power series in The ex¬ 
pression for Ug can be easily derived by means of the identity = x^G. After calculating 
6Ug = dllgllg^ one obtains the following explicit expression for B 2 , 


where 


B 2 


det 


sinhx 

X 


Ski + aPkOii bPkPi 
—bakOii 6ki — aakPi 


1 cosh X — 1 

2xsinhx 2x‘^ 2x‘^ 

1 1 cosh X — 1 

2a;sinhx 2x'^ 2x‘^ 


(37) 


(38) 


A general expression for the determinant with the structure of (^71) is given in Appendix 
A. The hnal result for B turns out to be remarkably simple 


B 


cosh a; 


X 


sinha; 


2V/+3 


(39) 


For A^j = 0 one immediately hnds that B = 1 and, for Nf = 1, the result is B = 
1 -|- ^( q ; i / 3 i -|- 0(2/32)• 

Instead of (^Tl) one could use an alternative parametrization dehned by 


U = VAV-\ 


(40) 


where A is a diagonal matrix with matrix elements Akk exp/^^, k = 1, • • •, Nf + 1 and 
AArj-i-2A7/+2 = exps. The matrix 1/ is a (compact) super-unitary matrix. The Berezinian 
from the transformation from the [/-variables to the V and A variables is given by 


nA:(e'^ — e*®fc)(e“^ — e“*®'=) 


(41) 


This parametrization is particularly useful if the integrals are calculated by means of a 
supersymmetric generalization of the Itzykson-Zuber integral |^, ^ However, in 
this case the Efetov-Wegner terms may be nonvanishing and must be carefully analyzed. 


5. The Microscopic Spectral Density in the Quenched Limit 


The valence quark mass dependence of the chiral condensate for Nf = 0 and ly = 0 was 
calculated in using the parametrization (^Tj). It is straightforward to generalize this 
computation to arbitrary topological charge. As a warm-up exercise for the calculation 
for arbitrary Nf, we hrst calculate S(m^) in a sector of topological charge z/ using the 
parametrization (^7|) for Nf = 0. 
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For Nf = 0 the parametrization (|^) yields a measure that is simply flat (see eq. (^)). 
In this genuinely quenched limit, the zero mode partition function (O) can be written as 

- J) = 

X 


d9 


271 


ds df3 da 


exp 


Sol/ Str 


TUy U 

0 Tlly — J 


(1 + |a/3) cos6* a(e*^ — e *)/2 
/3(e^ — e“*®)/2 (1 — ia/3) coshs 


The normalization of this partition function, Z’^^^^lrriy, my), for my ^ 0 follows by simply 
expanding the Grassmann variables. Using the Wronskian identity for the modified Bessel 
functions. 


K^{x)Iy+i{x) + Iy{x)Ky+i{x) = -, (43) 

X 

we And Z'^^^^irriy, my) = 1. This is in agreement with Wegner’s theorem for super-unitary 
invariant integrals |^. Moreover, in our context it has an immediate interpretation in 
terms of eq. (|l^, whose right hand side in this A^/ = 0 case simply reduces to a mass- 
independent constant (which conveniently can be chosen as unity, as done here). 

After differentiation with respect to J and a trivial integration over the Grassmann 
variables, the valence quark mass dependence of the chiral condensate is found to be 


S(m^) = 


de 


(cos e-cosh s) g^(ie-s) r cosh s(/i^ cos 0 + /i^ cosh s-1) 


(44) 


where Hy = myVEo. This integral is easily computed, it can be expressed in terms of 
Bessel functions. Using standard recursion relations and the Wronskian identity (|43|) the 
result we get is: 


T^irUy 






+ 


fly 


(45) 


The last term corresponds to the number of zero modes of the Dirac operator. It is 
remarkable that this term, which has such a simple interpretation in the framework of the 
original QGD partition function (where it comes from the explicit factor m'' in eq. (j|)), 
is reproduced by the supersymmetric chiral Lagrangian (which does not contain such an 
explicit factor m^). This phenomenon is general, and occurs in the framework of the usual 
chiral Lagrangian as well |^, 

The result (|45|) coincides exactly with the valence quark mass dependence of the chiral 
condensate obtained from Random Matrix Theory by integrating the microscopic spectral 
density according to mm- What is most important in the present context is that this 
relation can be inverted. Using eq. ([TT|) we immediately And 


u r 


Psid^) D Jvipjj) Ju+l{u)Jv—l{dJj) 


2 L 


(46) 
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This is the microscopic spectral density of the Dirac operator in QCD for iVj = 0 in a 
sector of arbitrary topological charge v. It agrees exactly with the result obtained earlier 
by means of chiral Random Matrix Theory |0. Here, it has been derived directly from 
the effective hnite-volume partition function of QCD in the microscopic scaling regime. 

The results derived in this section can be tested by means of lattice QCD simulations. 
Results obtained with staggered fermions for three colors convincingly show that the va¬ 
lence quark mass dependence of the chiral condensate [0 and the microscopic spectral 


density |^| are given by the above expressions for u = 0. Apparently, effects of the 
topological charge become only visible in simulations close to continuum limit. For com¬ 
pleteness we also mention that lattice QCD results for the microscopic spectral density 
for Nc = 2 with staggered fermions is given by similar universal expressions as well 


6. The Microscopic Spectral Density for QCD with Nf Flavors 


Let us now turn to the general case of Nf flavors and topological charge u. The 
(partially) supersymmetric effective partition function is given in eq. (|^). It will turn 
out that the measure of this partition function is normalized such that Z^^(A4) = Nf\ 
for J = 0 and zero sea quark masses. In general, it must satisfy the identity (^) up to a 
normalization factor. 

We compute the valence quark mass dependence of the chiral condensate by 


E{m,) = -djlogZ^^^{M) 


j=o 


(47) 


Using the parameterization ([27|) with U = UoUg, we can rewrite the partition function as 


= J d[U]SdeR(f/)exp|^Str[f/,Ad[/, + AdU-'f/-']|. (48) 


Separating out the fermion-fermion {FF) and the boson-boson {BB) blocks we have in 
the exponential 


Z’^f (M) = / d[[/]Sdet^(f/)exp SoU 




FF 

a 


M 


nUn + -MnU. 


FFtt-I 


U-^]- cosh S 


( 49 ) 


where the FF superscript stands for the upper left {Nf -|- 1) x [Nf + 1) block and BB 
represents the lower right element. The FF blocks of Ai and U are denoted by Ain and 
Un, respectively. We have also used the property that Ug and U~^, have the same FF 
blocks and BB blocks which is easily seen from the dehnition of Ug and its expansion in 
powers of the Grassmann elements. 
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Using that Sdet'^(U) = exp(—i/s)det'^(t/n) the s integral can be done immediately 
resulting in a modified Bessel function with argument = fij cosh, x where /ij is 

defined as /xj = (my — J)USo. The next step in evaluating the partition function is to 
expand the exponent in terms of the Grassmann variables in Ug^ = 1 + G^(cosha; — l)/x‘^ 
(see eq. (|35D ). Exploiting the rotational invariance of the scalar product x^, the exponent 
can be expressed in terms of the eigenvalues, A*, of the matrix 

A = EoViMnUn + U-^Mn)/2. (50) 

This allows us to make the following replacement in the integrand of (^) 

exp |Tr AUg^\ ^ exp Afc(l + —- aM 

I k ^ 

(51) 



Writing the exponent of the Grassmann variables as a product and once more using the 
symmetry of the Grassmann variables in the integral we find that an even more symmetric 
integrand is obtained by making the replacement 


exp 


f cosh X — 1 


x^ 


k > 


(1 — cosh 


z=o 

n 


= X! (1 — cosh 


X 


X 




n\ 


Aki ■ ■ ■ Aki 

l<ki<---<ki<n 

I (n-l)\ (dyY 


1=0 


n\ 


l\ 


y=0 


exp Tr log(l + yA). 

(52) 


The partition function can thus be written as 

yy P / '7^ \ 27T.~1“ 1 

= 2 J Y[d(3kdak cosh X 


n 

coshx) X! (1 — coshx)* 
/=o 


[n 


l)\ 


n\ 


fll(jCir 


(53) 


where 


^i(Mn) = / dUndetUn 


(dyY 


l\ 


det(l + yA) exp {TrA} 


(54) 


y=0 


In (q^ we have included our result for the Berezinian (|39|) which factorizes into an ordinary 
part an a Grassmannian part. 

Up to a factor (—I)""/!!, performing the final Grassmann integrals is equivalent to 
finding the coefficient of in the remaining integrand. Using the residue theorem the 
latter may be expressed as a contour integral around the pole at a; = 0 of the quotient of 
the integrand and We may then make the substitution z = sinhx to simplify the 

integrand resulting in 


ZUM) = 2 


dz (—l)’^n! 
27ri 


n 


I (n-l)\ 


1=0 


n\ 


ni(Mn). (55) 
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The term containing Ki, can be expanded using the product theorem for Bessel func¬ 
tions which can be written as 


+ z^) 


OO 1 

k=0 



(56) 


The other factor in the integrand can also be expanded in a series in z‘^. After multiplying 
the two out and extracting the coefficient of z‘^'^ we are left with the the hnal result for 
the partition function 


zuM)^2{-irY.j:T.i 


k=01=0 S=0 





(57) 


For 1 / = 0 this result can be further simplified by using the combinatorial identity 


5‘->- (0 (; 


U; (.p-i)¥ 


(58) 


for 0 < / < p and it gives zero otherwise. Here we take the right hand side of this 
expression to be one when p = / = 0. 

In general, the integrals over the unitary group give a fairly complicated expression. 
However, the calculation greatly simplifies for an arbitrary number of massless quarks 
and zero topological charge. Of course, the calculation is relatively simple for the case of 
Nf = 1, even with the quark masses included and arbitrary topological charge. 


7. Some Special Cases 


While the above formulas are quite involved in general, there are some important 
special cases in which they can be reduced to simple expressions in terms of elementary 
functions. Consider hrst the case of one physical quark of mass mi, i.e. Nf = 1. The 
expression ( |57[) then reduces to 

Z](r^=i(Ad) = ^iF;,(/ij)[z/(z/- 2)Ho - (i7 - -h H2] - ^AijA'i+^(/ij)[2i/Ho - Hi] 

+ (59) 

Using the Qi as computed in appendix B, we find that the normalization of the partition 
function is given by 


Zn,=i{M) 


J=0 




(60) 
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This result was to be expected. It coincides with the hnite volume partition function 
for Nf = 1 as derived in ^ (without valence quarks) and is just a special case of (p^Sl) . 
The valence quark mass dependence of the chiral condensate follows by differentiating 
with respect to the source term J. This results in 




u 


(/^u) T dj /_|_2 (p-t;)] T 

Pu 

„ .^;/(Pi;) Pu^fi/(Pi;)di/+1 (pi) pi/j/(pi)/j/-|_i(pj;) 

^- 




p? 


(61) 


This can easily be compared with results computed from Random Matrix Theory |^, ^ 
by integrating the microscopic spectral density obtained there according to eq. m- 
However, what is important here is not the precise form of this mass dependent chiral 
condensate, but again the fact that the relation ([l0|) can now be inverted to yield the 
microscopic spectral density itself. Using eq. (Q we hnd 

^ ^ It /'„.\2 7 (^.\ 7 7„,^l I Pl-^!^(pl)'^I/+l('*^) ~ 

Ps{u,mi) Jj/+i('u) h"! 7 / \ 2 I 2^ 

2L J U(pi) (u^ + pf) 


This result agrees exactly with the answer obtained from Random Matrix Theory P^ , PS 
One can of course continue by considering a larger number of massive flavors. But 
the computations rapidly get rather involved, and we have not succeeded in hnding a 
simple compact expression for the mass dependent chiral condensate for the general case 
of Nf massive flavors in a sector of arbitrary topological charge u. There is, however, 
one important case which can be computed rather easily: that of an arbitrary number of 
massless quarks in a sector of zero topological charge. In this case the calculation greatly 
simplihes since then = 0 for / > 3. This is readily seen from its dehnition in terms of 
the eigenvalues of A and the observation that for massless sea-quarks only two eigenvalues 
of A are nonzero. In this case the expression (|^ reduces to 


= 2 


KMOo+iY) I<„_i{pj)au 


n—1 


(63) 


where we have used the identity (see appendix B) 


(n - l)Ui + 2^2 = 0, 


(64) 


which is valid for the case of zero sea quark masses and topological charge equal to zero. 
From the explicit expressions in Appendix B for Hq smd Ui we obtain our hnal result for 
the partition function 


Zl=;{M) = {n-l)\ 


PjA^n.(pj) ZZTi ^ 


m,: 


m, 


m: 


n—2 


(65) 
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Next, from the Wronskian identity (|43|) we immediately find that Z(A4) = (n — 1)! = N/l 
for J = 0. The valence quark mass dependence of the chiral condensate follows from 
differentiation with respect to J, resulting in 


S(m^) 


— l^V 




( 66 ) 


This expression agrees with the result obtained previously from Random Matrix Theory 
[T8[| . Again our purpose here is just the opposite: we can now derive directly from field 
theory the microscopic spectral density of the Dirac operator from the discontinuity of 
S(m^) as given in eq. (^). The result is 


u r 




(u) = - 
K J 2 L 


JNfiuY — 


(67) 


This is the celebrated result obtained first using Random Matrix Theory Indeed, 
this taken together with the earlier results constitute an analytical proof that the smallest 
eigenvalues of the QCD Dirac operator are correlated according to a Random Matrix 
Theory whose form is dictated by the global symmetries of the QCD Dirac operator. 


8. Conclusions 


In the limit of vanishing light quark masses the infrared sector of the QCD partition 
function is dominated by the Goldstone modes associated with the assumed spontaneous 
breaking of chiral symmetry. However, the usual chiral Lagrangian does not allow us 
to access the Dirac spectrum. For that reason one has to extend the partition function 
with one valence quark and its superpartner |^. The chiral Lagrangian of this partition 
function is based on the super-group Gl{Nf + 1|1) x Gl{Nf + 1|1). In agreement with a 
supersymmetric generalization of the Vafa-Witten theorem and the maximum break¬ 
ing of chiral symmetry, the symmetry is broken to the diagonal subgroup Gl{Nf + 1|1). 
As dictated by the convergence of the integrals the integration manifold is restricted to 
a super-Riemannian symmetric submanifold. It is characterized by a symbiosis between 
compact and non-compact degrees of freedom. Remarkably, precisely this balance between 
compact and non-compact variables in the effective Lagrangian is what leads to the correct 
cut structure in the complex valence quark mass plane. The supersymmetric extension of 
the effective chiral Lagrangian has thus passed a highly non-trivial self-consistency test. 
Moreover, in the present context this super-extension of the chiral Lagrangian for QCD 
is not used to study artifacts of (partially) quenched numerical simulations, but rather to 
derive physical results in standard QCD with dynamical fermions. 

In a previous paper , three of us have shown that the microscopic spectral density 
can be obtained from a partition function with compact degrees of freedom only, i.e. by 
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replacing GI{1)/U{1) ^(1)- The integrals over the supergroups could be performed 

conveniently by means of a supersymmetric generalization of the Itzykson-Zuber integral. 
However, a similar approach applied to the direct calculation of the valence quark mass 
dependence of the chiral condensate only worked for the quenched approximation, i.e. for 
the G/(l|l) partition function, but failed in the case of a nonzero number of sea quarks. 
The general reason for such failure is well understood. The measure contains anomalous 
terms which have to be included. The appearance of these so-called Efetov-Wegner terms 
has its origin in the fact that the integration contour contains nilpotent terms, which after 
expansion in a power series, result in total derivatives. Typically, such terms contribute 
in the neighborhood of singularities. 

In this work we have followed a different approach. We have directly calculated the 
Gl{Nf + 1|1) partition function without relying on super-symmetric Itzykson-Zuber in¬ 
tegrals. Technically, this calculation was possible because we found a parametrization 
without anomalous contributions for the observables under consideration. 

In conclusion, we have shown analytically that the microscopic distribution of eigenval¬ 
ues of the QCD Dirac operator can be computed directly from a supersymmetric extension 
of the effective hnite-volume QCD partition function. The results agree exactly with the 
original computations which were based on chiral Random Matrix Theory. It is quite 
remarkable that what could appear as a forbiddingly difficult held-theory computation of 
the microscopic Dirac operator spectrum was first performed on the basis of universality 
arguments and Random Matrix Theory. What has now been proven is that these results 
are exact and indeed can be derived directly from held theory, without recourse to Ran¬ 
dom Matrix Theory. The underlying reason should be clear: Random Matrix Theories 
with the global symmetries of the QCD partition function can be reduced to the hnite 
volume partition function that has been studied in this paper. 
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Appendix A. Calculation of a determinant 

In this appendix we show that 


det ^ki + a(3kai h(3k(3i 
y —bakCXi 6ki — aakPi 


(1 — 2ax‘^ + (a^ — b‘^)x'^) 


( 68 ) 


Here, and (3k are Grassmann variables and a and b are scalar fnnctions. We use the 
notations that akCti represents an n x n matrix with entries akai and x"^ is dehned by 


X 


'y ^ (3k(^k- 


k=l 


a(3kai b(3k(3i 

—bakai —aak(3i 


The proof is as follows. If we introduce the matrix 

A = 

the determinant is calculated by the relation 

OO 

det(l + A) = ^ exp 


k=l 


(- 1 ) 


k+l 


k 


TrA^ 


(69) 


(70) 


(71) 


One can easily show that the powers of A have the same structure as the matrix A. 
Therefore, and 


_ f (3k^i bp f3kf3i 

y bp Q-kOli Qjp Otk(3i 

It is straightforward to derive a recursion relation for the coefficients Op and bp 



ttp+i = —ax^ttp — bx'^bp, (73) 

bp+i = —bx'^ttp — ax'^bp, (74) 

with oi = 1 and bi = b. The solution of these recursion relations is given by 

ap = ^{a + bYi-x^Y^-^'^+ ^{a-bY{-x^YP-^'^ 

bp = ^{a + bY{-xY^-^'>-^{a-bY{-x^)^-^^. (75) 

By resumming the power series into a logarithm and taking the trace of the matrices, one 
easily recovers the expression for the determinant. 
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Appendix B. Some integrals over Unitary groups 


In this appendix we calculate the integrals 


Qi(M)= / dUdetU" 


n 


det(l + yA) exp {TrA} , 


y=0 


where 


A = EoV(MU + U-^M)/2, 


(76) 


(77) 


and the integral is over the unitary group U{n). We evaluate these integrals for / = 0, 1 
and 2. These integrals are of the form 

.'.*■' (78) 


d[C/|/(C/, f/')det(f/)‘'exp + [/*)] ^ = (/((7,f/*) 


with M a diagonal matrix and U ranging over the group U{n). To evaluate these integrals 


we use the following result [47 


W{J,J^) = J d[U]det{Uy exp {Tt[JU + J^]} 


n—1 

n 2^A;! 

fc=i 


'det Jl’V detjj[z,^ ^Ij-i+^{zi)] 


det J 


A{zi 


(79) 


where Zi = 2a/\ with the Aj being the eigenvalues of the matrix and A{zf) is the 
Vandermonde determinant Ylk>i{^k ~ ^i) ■ From the invariance of the measure and the fact 
that an arbitrary complex matrix can be diagonalized by two unitary matrices resulting in 
a matrix with positive diagonal elements, it follows immediately that the integral factorizes 
into a function of the eigenvalues of and a ratio of determinants of and J. This 
result can be proven naturally [^ by means of a generalization of the Itzykson-Zuber 
integral to arbitrary complex matrices. We can now rewrite the original integral as 




J=Jt=M/2 


(80) 


Since W contains the eigenvalues A, we will need to change the derivatives with respect 
to the sources J and into derivatives on A by the chain rule 


d 


= E 


dXa d 


a=l ^^0 


(81) 


The derivative of A with respect to the sources can then be obtained from the identities 

n 

y a; = n(j<jy (82) 

p=l 
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for r = 1, • • - ,71. By applying derivatives with respect to the sources J and on both 
sides of the equation we get a system of equations which we can then solve for the required 
derivative of A. 

As an example we calculate the integrals corresponding to hli and ^2 as dehned in 
(ff^) which are used in the text. Our goal is to express them in terms of Oq =< 1 > and 
its derivatives with respect to the masses. We can then use the explicit expression for Oq 
for the two cases considered in the text to obtain an expression for these integrals. We 
consider the case with Nf = 1 sea quark with mass mi, one valence quark with mass m^ 
(n = 2) for topological charge u and the case with Nj massless flavors, one valence quark 
with mass m^ (n = Nf + 1) for zero topological charge (i/ = 0). In these cases, the values 
of rio are given by 


iA — ^-5-2- 

m„ — mf 


(83) 


and 


rio — 


j.r InA'^v) 

- Nf > 


Nf 

ruv 


(84) 


respectively. For simplicity we absorb the factors of Sold into the masses in this appendix. 
We emphasize that in this appendix the sea quark masses and the valence quark mass 
enter on the same footing. We nevertheless make this distinction since it allows us to 
make contact with the formulas in the main text. 

First we calculate 




=< TrA >= 


1 ” 
i 


d d \ 

94 ■^54/ 


hF( J, JA 


J=jt=M/2 ■ 


(85) 


One can check that the derivatives acting on the term (det det gives no contri¬ 
bution to the hnal result. Since this calculation only requires diagonal sources, it is easy 
to see that the derivatives acting on the remaining term gives 




E 


rrii 


d 

drrti 


Oq. 


( 86 ) 


For the two cases above, for Nf = 1 and one valence quark [n = 2) in the sector of 
topological charge z/ we obtain 


= 2y{my)y{mi) - 2Qo, (87) 

and for Nf massless flavors and one valence quark (n = Nf + 1) with zero topological 
charge (z/ = 0) we hnd 

= 2^fNf\^^A±l^^. (88) 

ruv ^ 
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The next integral, 


n2 = ^((TrA)2-TrA2), (89) 

includes off diagonal sources and thus becomes more complicated. The hrst term is similar 
to the previous example and can be written 


(Tryi)n= 


rrijm. 




_d _^ 

diTLi drrij 


■Q 


0- 


(90) 


Here again the term (det Jl/ det in W did not contribute. The second term is 

(TrH^) = 1 (t:t[2M^ + {MUf + {U ^. (91) 

We thus need to calculate < TiMUMU > and its conjugate. In order to evaluate these 
integrals we write 


TiiMUf) = ^ 






^_^ 

dJij dJji 






(92) 


and similarly for its conjugate. Now the term (det J^/ det does make a contribution. 
Summing the contribution for this case and its conjugate gives 2nv‘^VLQ. The remainder 
of these integrals can again be evaluated using the chain rule 


_a _ " d^K d " dXg dXb 

dJij dJji V SJijdJji d\a ^ 8Jij dJji dKdXb' 


(93) 


Applying the two derivatives to the identities (^) and then solving the resulting system 
of equations, we hnd (for i ^ j) 

dXa 


d Jijd Jji 


miTrij 

9 2 

— rUj 


(94) 


when i = a, i and j are reversed when j = a, and it is zero otherwise. Putting this all 
together we end up with 


d^n. 


0 


1 ” 

f!, = ' 4 


1 miTiij 

2 ^ - m] 


m. 


1 

4 

d^l, 

dm. 


n 

E 

i 

I dflo 
\^"dmi 

h 



m, 


>d^no 

dm"} 






(95) 


This can be further simplihed by using the identity 

mim, 


"W [to I'l J- 

2^ (2n-l)mj7^-hm,; 


' dmi 


dm"^ 


m'jQn 


+ 2E 


. m| - m] 


mi 


dflo 

dmi 


rrii 


d^o 

dmj 


= nuflr,. 


(96) 
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As a final result we obtain 


dQo lA d'^a 


0 mirrij 


- 2 (" - 1) S “‘8^ + 2 5 8^ - 5 


d^lr, diln 


2 2 I O O 

mj — 777.7 V orrii orrii 


( 97 ) 


For iVj = 1 with sea quark mass mi and valence quark mass the expression (|^) 
becomes for arbitrary u 


r22 = —+ Sz/hlo + ‘2myli+^{my)lu{my) , 


(98) 


while for Nf massless flavors and one valence quark we get for u = 0 


f^2 — 

2 •' 


(99) 


References 

[ 1 ] T. Banks and A. Casher, Nucl. Phys. B169 (1980) 103. 

[2] J.C. Osborn and J.J.M. Verbaarschot, Phys. Rev. Lett. 81 (1998) 268. 

[3] J. Gasser and H. Leutwyler, Phys. Lett. 188B(1987) 477; Nucl. Phys. B307 (1988) 
763. 

[4] H. Leutwyler and A. Smilga, Phys. Rev. D46 (1992) 5607. 

[5] E.V. Shuryak and J.J.M. Verbaarschot, Nucl. Phys. A560 (1993) 306. 

[6] M.A. Halasz and J.J.M. Verbaarschot, Phys. Rev. D52 (1995) 2563. 

[7] J. Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531; Phys. Lett. B329 (1994) 351; Nucl. 
Phys. B426[FS] (1994) 559. 

[8] G. Akemann, P. Damgaard, U. Magnea and S. Nishigaki, Nucl. Phys. B 487[FS] 
(1997) 721. 

[9] E. Brezin, S. Hikami and A. Zee, Nucl. Phys. B464 (1996) 411. 

[10] A.D. Jackson, M.K. Sener and J.J.M. Verbaarschot, Nucl. Phys. B479 (1996) 707. 

[11] T. Guhr and T. Wettig, Nucl. Phys. B506 (1997) 589. 

[12] A.D. Jackson, M.K. Sener and J.J.M. Verbaarschot, Nucl. Phys. B506 (1997) 612. 

[13] M.K. Sener and J.J.M. Verbaarschot, Phys. Rev. Lett. 81 (1998) 248. 


24 








[14] P.H. Damgaard, |hep-th/9807026 


[15] J.J.M. Verbaarschot, Lectures given at NATO Advanced Study Institute on Con¬ 
finement, Duality and Nonperturbative Aspects of QCD, Cambridge, 1997, Piep- 
th/971Uini. 


[16] J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, |hep-th/9806110| , accepted for 
pnblication in Nnc. Phys. B. 


[17] S. Chandrasekharan, Nucl. Phys. Proc. Snppl. 42 (1995) 475; S. Chandrasekharan 
and N. Christ, Nncl. Phys. Proc. Snppl. 42 (1996) 527; N. Christ, Nncl. Phys. B 
(Proc. Snppl.) 53 (1997) 253. 

[18] J.J.M. Verbaarschot, Phys. Lett. B368 (1996) 137. 

[19] J.J.M. Verbaarschot, in Nonperturbative Approaches to Quantum Chromodynamics, 
D. Diakonov, ed., Gatchina 1995. 


[20] P.H. Damgaard, Phys. Lett. B424 (1998) 322; G. Akemann and P.H. Damgaard, 
Phys. Lett. 432 (1998) 390; G. Akemann and P.H. Damgaard, Nucl. Phys. B528 
(1998) 411. 

[21] A. Morel, J. Physique 48 (1987) 1111. 

[22] C. Bernard and M. Golterman, Phys. Rev. D46 (1992) 853, Phys. Rev. D49 (1994) 
486; M. Golterman, Piep-lat/9411005| , Chiral Perturbation Theory and the quenched 
approximation of QCD. 


[23] B. Simons and A. Altland, |cond-mat/ 9811 134 


[24] M.R. Zirnbauer, J. Math. Phys. 38 (1997) 2007; J. Phys. A 29 (1996) 7113. 

[25] J. Verbaarschot and 1. Zahed, Phys. Rev. Lett. 70, 3852 (1993). 

[26] H. Leutwyler, Ann. Phys. 235 (1994) 165. 

[27] M. Zirnbauer, J. Math. Phys. 37 (1996) 4986. 

[28] K.B. Efetov, Supersymmetry in disorder and chaos, Gambridge University Press, 
(1997); K.B. Efetov, Adv. Phys. 32 (1983) 53. 

[29] F. Wegner, private communication, 1983. 

[30] M.R. Zirnbauer, Nucl. Phys. B 265 [FS15] (1986) 375. 

[31] T. Guhr, Nucl. Phys. A560 (1993) 223. 


25 















[32] M. J. Rothstein, Trans. Am. Math. Soc. 299, 387 (1987). 

[33] J. Alfaro and L.F. Urrntia, Piep-th/9810130 . 


[34] M. R. Zirnbaner and F. D. M. Haldane, Phys. Rev. B52 (1995) 8729. 

[35] Z. Plnhar and H.A. Wiedenmiiller, |cond-mat/9809276 . 

[36] T. Gnhr, J. Math. Phys. 32 (1991) 336. 


[37] J. Alfaro, R. Medina and L.F. Urrntia, J. Math. Phys. 36 (1995) 3085. 

[38] T. Gnhr and T. Wettig, J. Math. Phys. 37 (1996) 6395. 


[39] T. Gnhr, J. Math. Phys. 34 (1993) 2523; ibid. 34 (1993) 2541. 


[40] A.V. Andreev, B.D. Simons, and N. Tanignchi, Nncl. Phys. B432 [FS] (1994) 487. 

[41] P.H. Damgaard, Phys. Lett. B425 (1998) 151. 

[42] P.H. Damgaard, U.M. Heller and A. Krasnitz, [hep-lat/9810060| . M. Gockeler, H. 
Hehl, P.E.L. Rakow, A. Schafer and T. Wettig, Piep-lat/9811018| . 


[43] M.E. Berbenni-Bitsch, S. Meyer, A. Schafer, J.J.M. Verbaarschot and T. Wettig, 
Phys. Rev. Lett. 80 (1998) 1146; M.E. Berbenni-Bitsch, S. Meyer and T. Wettig, 
Phys. Rev. D58 (1998) 071502. 


[44] P.H. Damgaard and S.M. Nishigaki, Nncl. Phys. B518 (1998) 495. 

[45] T. Wilke, T. Gnhr and T. Wettig, Phys. Rev. D57 (1998) 6486; B. Seif, T. Wettig 
and T. Gnhr, |hep-th/9811044| . 

[46] G. Vafa and E. Witten, Nncl. Phys. B234 (1984) 173. 

[47] R. Brower, P. Rossi and G-I. Tan, Nncl.Phys. B190 [FS3] (1981) 699. 

[48] R. Brower and M. Nanenberg, Nncl.Phys. B180 [FS2] (1981) 221. 

[49] A.D. Jackson, M.K. Sener and J.J.M. Verbaarschot, Phys. Lett. B 387 (1996) 355. 


26 









